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Online Appendix I 



10 Calculation of the acceptance ratio for birth steps in bayou 

11 Consider a phylogeny with N tips and K shifts between selective regimes on a phylogeny with 

12 parameters &k = { a , °~ 2 , #o, #i> •••> s ii s /f, -^i, L K } where a and a 2 are the pull and diffusion 
is parameters of the OU process, 0j is the adaptive optimum for the ith regime (9 0 = root state), Sj is an 
14 index variable designating the branch on which the shift occurs, and Li is the location of the shift on 
is branch Sj. The probability of accepting a birth or a death move is 

R = min(l, LikelihoodRatio x PriorRatio x ProposalRatio x Jacobian) 

16 To move between parameter spaces with differing dimensionality, a branch is chosen at random, and if 

17 no shift is present, a birth step is proposed and a shift is added to the branch. If a shift is present, then 
is a death step proceeds and the shift is removed. Considering first the birth step, the adaptive optimum 

19 of the current branch, 9i, and descendant branches (up until a new rate is reached) are split into two 

20 rates, 9\ and 9' i+1 . This is done by first drawing a random number u <~ U(—v,v). In practice, we set v 

21 equal to 0.5, so that the proposal may be modified by a multiplicative tuning parameter S. The optima 

22 are then modified such that 

Q\ = 6 i — u5wi 

23 

=0i + u5w i+1 

24 where Wi is the weight assigned to the group inheriting optimum 9[ . All 6j where j > i are re-indexed 

25 as 9jj r \. 

26 Likelihood ratio. — Calculation of the likelihood for a multi-optima OU model has been outlined 

27 elsewhere (Hansen 1997, Butler and King 2004), and we therefore do not provide a derivation. 

28 However, in the most general terms the likelihood ratio is: 

p(Data\<d K+1 ) 
p(Data\e K ) 

29 Prior ratio. — The prior ratio is: 

prior(Q K+1 ) = P ( a ')p{or' 2 )p(K + l)p(P\K + 1) P (L'\K + l)p(9 Q )...p(9' l ),p(9' l+l )...p(9 K+l ) 
prior(e K ) ' p{a)p{a 2 )p{K)p{-t\K)p(Z\K)p{9 0 )...p{9 l )...p{9 K ) 
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30 Both p(s^\K + 1) and p(~t\K) are probabilities of a particular configuration of shifts among branches, 

31 which are simply 1 over the number of configurations K+l and K shifts can take among 2N — 2 

32 branches, respectively. Furthermore, only the parameters L' K+1 , Q\ and 6' i+1 are updated during a birth 

33 steps, with all other parameters remaining constant. As a result, the ratio can be simplified to: 
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35 Proposal ratio. — The proposal ratio is calculated to be: 
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where s* is a random variable that always take the value 1 and I is the proportional location of the 
shift on the proposed branch drawn from a uniform distribution [/(0, 1). The variable s* is included to 
maintain detailed balance between the forward and reverse moves. The denominator of the proposal 
ratio is the proposal probability for a birth from K to K + 1; which is the probability of proposing a 



birth move ( 



2N-2-K 



) times the probability of selecting that specific branch 



times the 



probability of selecting that specific value of u and I (the location of the shift on the branch) , given 
that u <~ U{—v, v) and I ~ U{1\, Z 2 ). The reverse is true for the numerator. Another way of thinking of 
it is as follows: Once a branch is switched from not having a shift to having a shift, it has an equal 
probability of being selected in the next move going in the reverse direction, and losing a shift. 



45 Jacobic 



The Jacobian is the matrix of partial derivatives: 



d(@ K +i) 



d(Q K ,s*,u,l) 

46 For applications in which the location of the shift can be placed in different parts of the branch, this 

47 variable would be uniformly distributed over the total length of the branch. 
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Acceptance ratio. — Putting it all together: 



A(Q K ,Q K+1 ) 



p(Data\e K+1 )p(K + 1)(K + l)p(e^)p(0' i+1 )(v 2 - v^jh - h)(wj + w i+1 )6 
p(Data\<d K )p(K){2N-2-K)p(6 t ) 



r' i+1 where r\ is the 



As is currently implemented, (vi,V2) = (—0.5,0.5) and Wi = r[ and w i+ i 
proportion of shifted total branch length inheriting 6[, r' i+1 is the proportion inheriting 0' i+1 . It follows 
that Wi + w i+ i = r[ + r' i+1 = 1. Furthermore, I2 — h = 1 since the proportional location of each shift is 
uniformly distributed along the branch [7(0, 1). Therefore, the the acceptance ratio becomes: 



A(e K ,e K+1 ) 
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Acceptance ratio for death steps 



Death steps are simply the inverse of the Jacobian and proposal ratio, with the subscripts changed as 



55 appropriate, to give: 



A(e K ,e K -i) = 



p(Data\0 K -i)p(K - 1)(2N - 1 - K)p(%) 
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Acceptance ratios for the rBM model of auteur 



58 The original version of auteur contained an error in the acceptance ratio for the addition or 

59 subtraction of a rate-shift [Eastman et al. 2011, eq. (3)]. We provide the corrected equation here, 
eo Consider a phylogeny with N tips and K number of rate shifts on a phylogcny with parameters 

ei 6 K = {y 0 , Si, s K ,r 0 , n, r K } where r$ is the ith rate, Sj is an index variable designating the 

62 location of the shift from rj_i — )■ r i7 and y 0 is the ancestral state. The probability of accepting a split 

63 or a merge move is 



64 To move between parameter spaces with differing dimensionality, a branch is chosen at random, and if 

es no shift is present, a birth step is proposed and a shift is added to the branch. If a shift is present, then 

66 a death step proceeds and the shift is removed. Considering first the birth step, the rate of the current 

e? branch, n , and descendant branches (up until a new rate is reached) are split into two rates, r[ and 

es r' i+1 . This is done by first drawing a random number u <~ U(—v,v). The rates are then modified such 

eo that 



R = mm(l, LikelihoodRatio x PriorRatio x ProposalRatio x Jacobian) 



r' i = ri — uSwi 
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r - +1 = n + u5w i+1 
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, where u>i is the weight assigned to the group inheriting rate r\ 



and 5 is a tuning parameter. 
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Likelihood ratio. — 
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p(Data\e K ) 
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Prior ratio. — The prior ratio is: 
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76 Proposal ratio. — The proposal ratio is calculated to be: 
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77 This follows as the denominator is the proposal probability for a birth from K to K+l; which is the 

78 probability of proposing a birth move ( 2J 2 jv-"^ ) times the probability of selecting that specific branch 

79 ( 2N-2-K ) t ne probability of selecting that specific value of u. The reverse is true for the 

so numerator. Another way of thinking of it is as follows: Once a branch is switched from having a shift 

si to not having a shift, it has an equal probability of being selected in the next move of going in the 

82 reverse direction, and losing a shift. Thus, there is no correction necessary for the different proposal 

as probabilities. 

84 Jacobian. — The Jacobian is the matrix of partial derivatives: 

d(e K +i) 



d(@ K ,s*,u) 

as where s* is a random variable that always take a of value 1. For applications in which the location of 

86 the shift can be placed in different parts of the branch, this variable would be uniformly distributed 

87 over the total length of the branch. 
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on ^ccepiance roiio. — Putting it all together: 
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91 If, as is currently implemented, (^1,^2) = (— ^nj, r i+:L n i+ i) and tUj = and w i+ i = n i+ i where nj is 

92 the number of branches inheriting r-, then the acceptance ratio becomes: 
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Online Appendix III 

94 Based on the results our preliminary simulation studies, we found that the number of shifts 

95 seemed particularly sensitive to specification of the prior. We conducted several analyses to establish 

96 that the method is capable of using information in the data to adequately estimate the number of 

97 shifts. In addition, we conducted a prior sensitivity analysis to evaluate the effects of prior specification 

98 on biological inference. 

99 Methods 

100 Balanced tree test case. — We tested whether our method performed as expected when confronted with 

101 a simple test case in which the number of optima should be easily inferred. We simulated a balanced 

102 phylogeny with 64 tips with equal branch lengths and a total tree height scaled to one. We set a at a 

103 very high value of 10 and a 2 = 0.1. We placed four shifts on the phylogeny resulting in four clades of 

104 16 tips of equal age with optima of 9 = {—2, —1, 1, 2} and a root state of 0. Under these conditions, 

los very little covariance exists between tip states and the distribution resembles an i.i.d mixture of normal 

loe distributions, for which our method should be expected to find the correct number of shifts. We ran 

107 bayou with a strong prior on only a few shifts (conditional Poisson with A = 1, K max = 32) using 4 

los indcpcdcnt chains for 500,000 generations. We combined the resulting chains after discarding the 

ioo burn-in and analyzed the resulting posterior of shift numbers. 

no Bayes Factor comparisons using fixed numbers of regimes. — Estimation of the number of shifts is 

in complicated by the very large number of models that are accessible by the reversible-jump algorithm 

112 which may have relatively equivalent likelihoods. To simplify our exploration of this space and 

n3 establish whether the number of shifts was reasonably estimated, we sought to test the model 

n4 performance across a single slice of marginal likelihood surface. To do so, we simulated a phylogeny 

us with 64 tips. We arbitrarily selected 5 clades with at least 5 tips and assigned each to a unique regime. 

lie We then simulated the data with parameters a = 3, a 2 = 1 and 8 0 = 0. The remaining optima 

H7 {6i, 9 5 } were drawn from a normal distribution with mean = 0 and sd = 5. We then set up priors in 

us which the location of the regimes were fixed, thus eliminating the reversible-jump aspect of the 

no algorithm and simplifying analysis (while also speeding convergence). We fit both simplified and more 

120 complex versions of the true model (with 5 shifts) by including only a subset of the true shifts or by 

121 adding additional shifts, respectively. We thus fit models with K — 0, 1, 2, 3, 4, 5 (i.e. the true model), 6 

122 or 10 shifts. These models are therefore nested, with each model being a more general model than the 

123 previous one. For each model, we ran bayou for 200,000 generations and estimated an initial Markov 
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124 chain. We then used these Markov chains to initialize the stepping stone algorithm to estimate the 

125 marginal likelihood, as described in the main text. We estimated power posteriors in 5 steps along the 

126 sequence from 0 to 1, with each chain run for 200,000 generations. We then computed Bayes Factors 

127 comparing the true model with 5 shifts to the reduced (K < 5) or overparameterized {K > 5) models. 

128 Prior sensitivity analysis. — We simulated 20 phylogenies with 128 terminal taxa under a pure-birth 

129 model with a birth rate of 0.1 as in the main text. Trees were scaled to unit height. For each 

130 phylogeny, four different regime configurations were painted onto the branches corresponding to K= 1, 

131 5, 10, or 20 shifts. The magnitude of each shift is drawn from a uniform distribution. The minimum 

132 and maximum values of the uniform distribution corresponded ±8y/ c 2 / (2a). Given that a 1 and a 

133 were both set to 3 (as in the main text), the shifts ranged from ±5.66 in magnitude. Shifts were spaced 

134 equally through time according to the number of shifts in the simulation. If only one shift was present, 

135 it was placed 0.7 tree units before present. Otherwise, K shifts were spaced evenly across the interval 

136 from 0.9 to 0.1 tree units before present and were placed on randomly selected branches present at a 

137 given time point. This was done to aid in examination of the effect of age and magnitude on the 

138 estimation of shift number and location. 

139 Simulated datasets were then analyzed using one of three prior distributions for the total 
no number of shifts, either a conditional (i.e. truncated) Poisson distribution, a discrete uniform 

i4i distribution or a negative binomial distribution. A total of 13 different parameter sets were chosen, 5 
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Figure 1: Prior densities used in the prior sensitivity study. Three different distributions were used, the 
conditional Poisson (top), the discrete uniform (middle) and the negative binomial (bottom). 
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142 for the conditional Poisson distribution (A = {1, 5, 10, 15, 20}, K max = 64), 2 discrete uniform 

143 distributions ({min = 0, max = 64}, {min = 1, max = 32}) and 6 negative binomial 

144 ({size = l,prob = 0.1}, {5, 0.357}, {50, 0.847}, {5, 0.833}, {5, 0.5}, {5, 0.2}) which correspond to means 

145 of {9, 9, 9, 1 ,5, 20} and variances of {90, 25.2, 10.6, 1.2, 10, 100} respectively. Priors for other 

146 parameters were the same as in the base simulations of the main text: a ~ LogNormal(ln fi = 0.25, In 
i4 T a = 1.5), a 2 ~ LogNormal(ln fj, = 0, In a = 5), 6 ~ Normal (/x = 0, a = 3). Each MCMC was 

148 estimated using bayou using two independent chains run for 500,000 generations, with the first 0.3 of 

149 each chain being discarded as burn-in. Only runs in which the Gelman and Rubin's R-value for all 

150 parameters reached below 1.1 were included in subsequent analyses. 

151 Results 

152 Balanced tree test case. — All chains eventually converged on at a minimum of three shifts, which is the 

153 simplest model that allows each of the four clades to occupy a unique optimum. The placement of 

154 these shifts varied across chains, and runs did not converge for a owing to the very large range of 

155 values that a can take (Supplementary Figure [2J). Furthermore, chains often found alternative local 

156 optima for the location of shifts, and mixing was slow between these alternative configurations of shifts, 

157 consistent with our previous observation that steep likelihood peaks result in poor convergence and 

158 inadequate mixing. Nevertheless, there was very strong support for high a, distinct optima in each 

159 clade and a minimum of four phenotypic optima, consistent with expectations. 

160 Bayes factor comparisons using fixed numbers of regimes. — The best fitting model was the true model, 




Figure 2: Estimation of shifts on a balanced phylogeny under ideal conditions. (A) Posterior support is 
indicated by the size of the circles across four combined chains. (B) A phonogram showing the distribution 
of traits across the phylogeny. The solid red line is the posterior density for the phenotypic optima (9) 
from the combined MCMC chains. (C) Posterior of the number of shifts compared to the prior. All runs 
converged on a solution with at least 3 shifts, giving each clade a unique optimum (including the root). 
Note that these runs did not converge for parameters a and a 2 . 
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Figure 3: Estimation of Bayes Factors across an axis of increasing complexity relative to a simulated 
model. (A) A phenogram showing the true location and number of shifts used to generate the data 
(B) The maximally complex model tested, with 5 additional shifts added at arbirtrary locations (C) 
Marginal likelihood surface across models with increasing numbers of shifts. True shifts were added 
randomly one at a time until the true model was reached (K = 5, dotted line), after which arbitrary 
shifts were added to increase the complexity of the model. (D) 2 In BF indicating the support for the 
true model against all other models. The true model had strong support against simplified as well as 
more complex models. Support was stronger for the true model when compared to simplified models 
than it was for more complex models. 

lei with all 2lnBF > 6, indicating strong support (Kass and Raftery, 1995). This indicates that there is 

162 often substantial information to inform the number of shifts under our model. The true model was 

163 more strongly favored over the simpler models than it was over more complex (Supplementary Figure 

164 [3]). This asymmetry suggests that the method has more difficulty distinguishing overparameterized 

165 models from the true model relative to underparameterized models. 

lee Prior sensitivity analysis. — Our prior sensitivity analysis demonstrates that the estimated number of 

167 shifts is sensitive to the prior distribution. Of the three distributions we tested, the conditional Poisson 

168 distribution influenced the posterior distribution most heavily (Supplementary Figure [4| . The number 
leg of shifts tended to be underestimated in simulated datasets when K = 10 (except for the conditional 

170 Poisson with A = 20), and all priors severely underestimated the number of shifts when K = 20. In 

171 fact, across priors, most datasets with K = 20 resulted in as many shifts being identified as when 

172 K — 1 (Supplementary Figure [4]). 

173 To explore the reasons for this behavior, we examined model fits for each of the four numbers 

174 of simulated shifts. We find that when the number of shifts is underestimated for the 20 shift datasets, 
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175 it is because a large number of runs converged on a values that are very low (Supplementary Figure [5]). 

176 Simulations in which the estimated a value approached the true value of 3 came much closer to 

177 correctly estimating the number of shifts. We conclude that when the number of shifts is high, it is 

178 difficult for the method to distinguish between shifts resulting from optima shifts, and within-regime 

179 divergence. Thus, trees with many shifts begin to resemble BM-likc models. Under these conditions, 

180 the method is likely to simply return the prior on the number of shifts. This is because when a is low, 
lsi any number of shifts can occur without substantially affecting the likelihood of the model. 

182 While the number of shifts may be sensitive to the prior, we were interested in whether the 

183 overall number would dramatically affect inference of particular shifts on branches. As in the main text, 

184 older and higher magnitude shifts are well estimated across all K — 1, 5 or 10 models (Supplementary 
las Figures [6] and [7| . Branches containing shifts that return low posterior support tend to have shifts that 

186 are recent and of low magnitude. This helps explain why the number of shifts tends toward being 

187 underestimated when the number of shifts is large: recent shifts become more common (due to most 
las branches being near the tips) and are unlikely to be detected by the method. We show that the branch 

189 posterior probabilities of having a shift are relatively insensitive to the prior (Supplementary Figures [6] 

190 and [7] and correlations between branch posterior probabilities across prior sets are extremely high 

191 (Supplementary Figure [8]) with the exception of 20-shift models, which have low correlations across 

192 branch lengths. These low correlations were primarily the result of few 20-shift datasets having any 

193 highly supported shifts, resulting in overall very low posterior probabilities across all branches. 

194 We illustrate these effects using several exemplar simulations. When the shift is of high 

195 magnitude, all prior sets result in high posterior probabilities for the true shift while all other posterior 

196 probabilities remain low (Supplementary Figure [9b. When the magnitude of the shift is small, as in 

197 Supplementary Figure [10[ under all priors the method unsurprisingly fails to detect a high posterior 

198 probability for the true shift. We note that in both cases, regardless of the estimate of the total number 

199 of shifts, branches without shifts have low posterior support. For example, for in both the high and low 

200 magnitude case, the prior cpois ~ (A = 20) returns an estimated number of shifts of around 10. Yet no 

201 individual branches show dramatically elevated posterior probabilities unless a true shift is present. 

202 When the true number of shifts is high, many of these shifts may be missed (Supplementary 



Figure 11 & 12). Those that are missed tend to be of low magnitude and recent age (Supplementary 



Figure 111. In many cases, only by imposing a very strong prior for a high number of shifts (cpois 
~ A = 20) were we able to obtain support for shifts in optima across the phylogeny and approach the 



true parameter values of the simulated model (Supplementary Figure 11). In some cases, the entire 
model collapsed and became " BM-like" , with no shifts identified and estimates of a effectively equal to 



0 (Supplementary Figure 12). Thus, instead of resulting in false positives, the method conservatively 
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209 returns simplified models, even when the prior on the number of shifts is high. 

210 Discussion 

211 Given the results of the prior sensitivity analyses, we conclude that our method can reliably estimate 

212 the location and magnitude of shifts in adaptive optima as long as the number of shifts is not large, 

213 and the prior allows sufficient complexity to fit the data. Inference should focus on the branch posterior 

214 probabilities, not the total number of shifts returned by the model, which can be sensitive to prior 

215 specification. By contrast, the branch posterior probabilities for particular branches with shifts tend to 

216 have elevated support regardless of the prior, while branches without shifts had low support across 

217 methods even when the number of estimated shifts was high. Furthermore, the results demonstrate 

218 placing a strong prior that undcrparameterizes the model results in worse model-fits than when the 

219 model is overparameterized by placing a strong prior with large number of shifts. We find that if the 

220 true number of shifts is much larger than the specified prior expectation, the model will often collapse 

221 into a BM-like model in which the estimated shifts are explained not by shifts in optima, but by 

222 incorporating this divergence into the parameter a 2 . By contrast, specifying a prior with a large 

223 number of shifts is relatively benign and will only very rarely result in specific branches having high 

224 posterior probability when no shift is present. Instead, all branches will have only slightly elevated 

225 posterior probabilities, shifts will be reconstructed to be of low magnitude, and estimation of other 

226 parameter values such as a 2 and a will be only slightly affected. We therefore recommend against using 

227 conservative priors with a low expected number of shifts, and obtained the most reliable inference 

228 across simulation conditions when using a conditional Poisson prior with a high A. 

229 Identifying complex models is difficult, and in many cases in our prior sensitivity study, it may 

230 have been possible that given enough generations, our MCMC's could have converged upon the true 

231 model. Certainly, some recent, low magnitude shifts will never be identified as having high posterior 

232 support, and in this sense, inference on the true number of shifts is a poorly formed question. Rather, 

233 our method should be used to infer the location of shifts of sufficient magnitude and age that have 

234 affected the distribution of trait data across the phylogeny. By default, bayou draws starting values at 

235 random from the supplied prior distribution, and these can have very low likelihoods. Consequently, 

236 runs often quickly remove most shifts from the phylogeny and slowly add in complexity through the 

237 course of the MCMC. This often leads to the MCMC exploring parameter space in regions with very 

238 low a values, which in turn results in very little gain in likelihood support even when shifts are 

239 correctly identified by a birth proposal. An alternative strategy is to supply bayou with an 

240 overparameterized model as the starting point (e.g. each tip has a unique shift with an optimum near 
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Figure 4: Distribution of posterior mean estimates of the number of shifts from phylogenies with K = 1 
(top left), K = 5 (top right), K = 10 (bottom left) and K = 20 (bottom right) true shifts. Note, this is 
the distribution of the mean estimated number of shifts across 20 independent simulations, rather than 
the posterior distribution itself. 



241 its trait value). This ensures that complex models remain accessible to the MCMC, while still allowing 

242 the method to converge on BM-like or single optimum OU models, which can be found easily regardless 

243 of the starting point. However, when the inferred models are complex with high a values and a large 

244 number of shifts, several chains should be run from different starting points to ensure that the MCMC 

245 does not getting stuck on local optima. 
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a a 

Figure 5: Distribution of estimates for a and the number of shifts K across the prior sensitivity analysis. 
Each point represents the estimate of a and K for a single simulation, with the color indicating the 
prior used. Dotted lines indicate the true parameter values, solid line is the density of a values across 
all simulations and priors. Mis-estimation of a appears to drive mis- estimation of the number of shifts. 
In particular, note that for when K = 20 in the true model, most models collapse to a values near 0, 
indicating BM-like evolution. For these cases, only those models with a high prior on the number of 
shifts (cpois ~ A = 10 and cpois ~ A = 20) approach the true model. 
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Figure 6: Posterior probabilities for individual branches containing a shift. Each panel contains all 
branches for 20 trees with the indicated number of true shifts simulated. Each line indicates the posterior 
support for a shift occuring on a single branch estimated across the 13 different MCMC runs using each 
of the different prior distributions. Shifts are colored by their magnitude multiplied by their age before 
present. Higher values indicate older, higher magnitude shifts. Gray points indicate the average posterior 
value of branches without shifts. Note that older, higher magnitude shifts have high support across all 
priors when the true model contains only a few shifts. However, when 20 shifts are present, most branches 
contain low posterior support, even when high in magnitude and age. 
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Figure 7: The same data as in figure |6] but averaged over shifts of equal age. Gray lines are the average 
across branches without shifts. Again, notice that there is lower posterior support for all branches when 
the number of shifts is large. 
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Figure 8: Correlation in branch posterior probabilities across priors when the true model contains 1 shift 
(A), 5 shifts (B), 10 shifts (C) or 20 shifts (D). Correlations are high regardless of the prior used until the 
model collapses into BM-like evolution (see Supplementary Figure [5]). This occurs only when the overall 
variation in branch posterior probabilities across branches is small (no shifts are highly supported). 
Furthermore, some simulations recovered high correlations even for the case of 20 shifts (D), but are not 
apparent in the aggregate data. 
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Figure 9: Example simulation and estimated parameters for a phylogeny with a single shift. (A) 
Phenogram showing the distribution of traits and location of true shifts. (B) Estimated a and K 
across prior distributions (true values indicated by dotted lines). (C) Relationship between shift age and 
magnitude and the branch posterior probabilities (D) Density of branch posterior probabilities. Colors 
for plots B-D indicate prior distribution. Note that while some prior distributions over estimate K, the 
posterior probabilities of individual branches and the estimation of a are only slightly affected. Only the 
branch containing the true shift shows a high posterior probability, with little variation across priors. 
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Figure 10: Example simulation and estimated parameters for a phylogeny with a single shift in which the 
shift is of low magnitude. (A) Phenogram showing the distribution of traits and location of true shifts. 
(B) Estimated a and K across prior distributions (true values indicated by dotted lines). (C) Relationship 
between shift age and magnitude and the branch posterior probabilities (D) Density of branch posterior 
probabilities. Colors for plots B-D indicate prior distribution. Again, while K is overestimated for some 
priors, posterior support for individual branches is low across priors. 
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Figure 11: Example simulation and estimated parameters for a phylogeny with 20 shifts. (A) Phonogram 
showing the distribution of traits and location of true shifts. (B) Estimated a and K across prior dis- 
tributions (true values indicated by dotted lines). (C) Relationship between shift age and magnitude 
and the branch posterior probabilities (D) Density of branch posterior probabilities. Colors for plots 
B-D indicate prior distribution. Most prior distributions failed to recover any shifts. The best-model at 
recovering true shifts was the cpois ~ A = 20 followed by cpois ~ A = 10. All other models collapsed to 
BM-like models with low a values (some priors are not shown because runs failed to converge) . Further- 
more, only high magnitude shifts of sufficient age were strongly supported while most low magnitude, 
recent shifts were missed (C) and the number of shifts was underestimated (B). 



20 




Figure 12: Example simulation and estimated parameters for a phylogeny with 20 shifts. (A) Phonogram 
showing the distribution of traits and location of true shifts. (B) Estimated a and K across prior 
distributions (true values indicated by dotted lines). (C) Relationship between shift age and magnitude 
and the branch posterior probabilities (D) Density of branch posterior probabilities. Colors for plots B-D 
indicate prior distribution. All prior distributions failed to recover the true model or strongly support 
any branches (C & D), even when the total number of estimated shifts was high (B). 
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